Light from Cascading Partons in Relativistic Heavy-Ion Collisions 
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We calculate the production of high energy photons from Compton and annihilation processes as 
well as fragmentation off quarks in the Parton Cascade Model. The multiple scattering of partons 
is seen to lead to a substantial production of high energy photons, which rises further when parton 
multiplication due to final state radiation is included. The photon yield is found to be directly 
proportional to the number of hard collisions and thus provides valuable information on the pre- 
equilibrium reaction dynamics. 



High-energy collisions of heavy nuclei are studied with 
the goal to explore the properties of matter at energy den- 
sities many times that of normal nuclear matter. There 
are strong theoretical reasons to believe that matter un- 
der these conditions, when thermally equilibrated, forms 
a Quark-Gluon Plasma (QGP), a state of matter in which 
quarks and gluons are not confined to hadrons. The dis- 
covery of this state of matter and the study of its proper- 
ties would constitute a major advance in our understand- 
ing of quantum chromodynamics (QCD). 

The parton cascade model (PCM) |lj was proposed to 
provide a detailed description of the temporal evolution 
of nuclear collisions at high energy, from the onset of hard 
interactions among the partons of the colliding nuclei up 
to the moment of hadronization. The PCM is based on 
a relativistic Boltzmann equation for the time evolution 
of the parton density in phase space due to perturbative 
QCD interactions including scattering and radiation in 
the leading logarithmic approximation. 

During the epoch where the partonic picture of sequen- 
tial perturbative interactions among QCD quanta may 
be valid, important many-body effects, such as parton 
shadowing, color screening, and Landau-Pomeranchuk- 
Migdal (LPM) suppression of gluon radiation, complicate 
the single-particle dynamics of parton transport. Obvi- 
ously, the PCM involves the application of perturbative 
QCD to a hitherto untested domain. This poses the ques- 
tion - how can one test the correctness of this approach? 
This is clearly not possible by studying the hadronic fi- 
nal state, since - if the dense matter equilibrates - all in- 
formation about its dynamics before equilibrium will be 
lost. One may look for remnants of the pre-equilibrium 
dynamics in the form of very energetic partons, which 
did not lose a large fraction of their initial momentum 
due to interaction in the dense medium and produce jets. 
Even better suited to this purpose are weakly interacting 
probes of the pre-equilibrium dynamics, such as quanta 
that interact only clcctromagnctically: photons and lep- 
tons. 

Photons, which are produced from Compton (qg — > 
qj), annihilation (qq — » 57), and bremsstrahlung (q* — > 
97) processes (see Fig. EJ, thus assume considerable im- 



portance. These are analogous to the the processes gov- 
erning the energy loss of energetic partons, where gluons 
are emitted instead of photons. Photons, however, once 
emitted, almost always escape the system, and their en- 
ergy and momentum will remain unaffected Q as they 
will only rarely interact with the surrounding medium. 
For this reason they are excellent and reliable probes of 
the evolution of the partonic cascade in nuclear collisions. 
In fact, one can argue that the photons confirm the pres- 
ence of parton cascading processes after the initial pri- 
mary parton-parton collisions. 

In this letter we address several questions relating to 
the value of photons as probes of the dense matter cre- 
ated in a nuclear collision: Arc photons sensitive to the 
pre-equilibrium partonic reaction dynamics, i.e. parton 
re-scattering? What effect do parton fragmentation pro- 
cesses have on photon production? Is there a direct 
connection between the photon yield and the number 
of hard perturbative parton-parton collisions during the 
reaction? In view of our specific focus, we do not con- 
sider the production of photons from a late-stage, ther- 
mal mixed phase or hadronic phase. These contributions 
are quite well understood 3] , and can easily be added to 
our results. 

The fundamental assumption underlying the PCM is 
that the state of the dense partonic system can be 
characterized by a set of one-body distribution func- 
tions Fi(x IJ ',p a ), where i denotes the flavor index (i = 
g, u, u, d,d, . . .) and x M ,p Q are coordinates in the eight- 
dimensional phase space. The partons are assumed to 
be on their mass shell, except before the first scatter- 
ing. In our numerical implementation, the GRV-HO 
parametrization 4] is used, and the parton distribu- 
tion functions are sampled at an initialization scale Qo 



0= 



(p™ m ) 2 ; see later) to create a discrete set of particles. 



Partons generally propagate on-shell and along straight- 
line trajectories between interactions. Before their first 
collision, partons may have a space-like four-momentum, 
especially if they are assigned an "intrinsic" transverse 
momentum. 

The time-evolution of the parton distribution is gov- 
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FIG. 1: Photon production processes included in our calcula- 
tion. 
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erned by a relativistic Boltzmann equation: 

P^F^P)=Ci[F] (1) 

where the collision term Ci is a nonlinear functional of 
the phase-space distribution function. The calculations 
discussed below include all lowest-order QCD scattering 
processes between massless quarks and gluons |5| , as well 
as all (2 — > 2) processes involving the emission of pho- 
tons {qg -> qj, qg -> 37, qq -> gj, qq -> 77). A low 
momentum-transfer cut-off p™ m is needed to regularize 
the infrared divergence of the perturbative parton-parton 
cross sections. A more detailed description of our imple- 
mentation is in preparation Q. 

We account for the final state radiation Q following 
a hard interaction in the parton shower approach. In 
the leading logarithmic approximation, this picture cor- 
responds to a sequence of nearly collinear 1 — > 2 branch- 
ings: a — > be. Here a is called the mother parton, and 
b and c are called daughters. Each daughter can branch 
again, generating a structure with a tree-like topology. 
We include all such branchings allowed by the strong and 
electromagnetic interactions. The probability for a par- 
ton to branch is given in terms of the variables Q 2 and 
z. Q 2 is the momentum scale of the branching and z 
describes the distribution of the energy of the mother 
parton a among the daughters b and c, such that b takes 
the fraction z and c the remaining fraction 1 — z. The 
differential probability to branch is: 

b,c 

where the sum runs over all allowed branchings. The 
a a bc is a em for branchings involving emission of a pho- 
ton and a s for the QCD branchings. The splitting ker- 
nels P a ~*bc are given in 0\. The collinear singularities 
in the showers are regulated by terminating the branch- 
ings when the virtuality of the (time-like) partons drops 
to /!q, which we take as 1 GeV. In principle, one could 



FIG. 2: Transverse momentum distribution of photons from 
central collision of gold nuclei at RHIC (see text for an expla- 
nation of the symbols). 



take a smaller value for the cut-off fi for a quark frag- 
menting into a photon 9], but we have not done so as 
we are only interested in high energy photons here. The 
soft-gluon interference is included as in 0, namely by 
selecting the angular ordering of the emitted gluons. An 
essential difference between emission of a photon and a 
parton in these processes is that the parton encounters 
further interactions and contributes to the build-up of 
the cascade, while the photon almost always (in our ap- 
proximation, always) leaves the system along with the 
information about the interaction. 

Some of these aspects were discussed by authors of 
Ref. 0| and explored within framework of the PCM im- 
plementation VNI 12]. Unfortunately, several compo- 
nents of their numerical implementation were incorrect. 
Our present work is based on a thoroughly revised and 
extensively tested implementation of the parton cascade 
model by the present authors, called VNI/BMS 0. 

As a first step, we discuss the results obtained when 
the parton cascade is restricted to two-body scattering. 
We are considering central collisions (b — 0) of Au nu- 
clei at full RHIC energy (y/s = 200 GeV per nucleon 
pair). We have confirmed the accuracy of our code by 
comparing the analytical pQCD prediction for prompt 
photon production, due to Compton and annihilation 
processes, with our results, when we use an eikonal ap- 
proximation in the cascade. Our results are shown in 
Fig. [21 where we compare the contributions from inter- 
actions involving at least one primary parton (triangles) 
with the result obtained with full binary cascading (di- 
amonds). Note that contribution from primary partons 
does not reflect the momentum broadening (Cronin ef- 
fect) due to multiple soft scattering. The figure shows 
that most of the photons between 2 and 4 GeV have 
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FIG. 3: Upper frame: photon transverse momentum spectra 
for contributions from collisions and fragmentations (trian- 
gles) and from collisions alone (squares), when parton frag- 
mentations are included in the calculations. Results for the 
calculations with only the binary collisions are also given (di- 
amonds). Lower frame : The multiplicity density of photons 
having pr > 2 GeV. 



their origin in the multiple semi-hard scattering of par- 
tons. This lends support to a recent prediction that there 
would be a substantial production of photons from the 
passage of high energy quarks through the QGP which 
may be produced in such secondary collisions jl-'if . We 
find the average squared transverse momentum of the in- 
coming rescattered partons at the photon vertex to be 



2.1 (GeV/c) 2 , confirming previous phenomcno- 



logical analyses of nuclear broadening effects in prompt 
photon production (l4|. 

We next allow for fragmentation processes after binary 
interactions as discussed above. The lower set of points 
in the upper frame of Fig. |21 (squares) includes only QCD 
induced fragmentations of partons. A comparison with 
Fig. |21 (diamonds) shows that this already leads to an 



increase in the photon emission by about a factor 4. The 
reason for this surprisingly dramatic effect is that the 
fragmentation leads to a large multiplication of partons 
and thus to an enhanced photon yield to the increase in 
the number of multiple scatterings. In comparison, the 
increase in the photon yield due to the fragmentation 
of photons off the scattered quarks (triangles) is small, 
except below 1 GeV, filling up the artificial gap caused 
by the cut-off p™ m in the binary collisions. The strong 
enhancement in the low pt region is not surprising, since 
radiative processes generally favor low-energy final states. 

The lower frame of Fig. shows the rapidity distribu- 
tion of photons with pr > 2 GeV for the same two cases 
as are discussed in the upper frame, but includes also 
the distribution obtained when only binary interactions 
are allowed (diamonds, corresponding to the uppermost 
curve in Fig.[5J. The enhancement of photon production 
due to the inclusion of fragmentation processes occurs 
mainly around mid-rapidity, where a dense partonic sys- 
tem is created and most of the fragmentation processes 
contribute. We chose the cut-off px > 2p^, because the 
photon yield below this value is expected to receive a 
large contribution from hadronic reactions, which do not 
interest us here. 

The strength of the source of single photons is often 
discussed in terms of the ratio (7/71-°), as most of the 
back ground to them comes from the decay of ir° . We can 
not get a precise value for this unless we supplement our 
description with a model for the hadronization. A reason- 
able alternative is provided by the ratio of (7/partons) 
as the multiplicity of pions will be proportional to the 
multiplicity of partons. We find that this quantity is of 
the order 0.04% when only multiple scatterings among 
partons is considered and rises to 0.15% when fragmen- 
tations are allowed and which considerably enhance the 
multiple scatterings. 

The results shown so far do not account for the pos- 
sible LPM suppression of radiation from the scattered 
partons. An accurate implementation of this effect in a 
semi-classic treatment is difficult. To explore the pos- 
sible magnitude of this effect, we consider here a treat- 
ment in which the scattered parton, having energy E a 
and time- like virtuality m 2 , is forbidden to fragment if 
it re-scatters before the formation time of the radiated 
parton- E a jm? a has expired. This ansatz should provide 
us with an upper limit of the suppression as it does not 
account for the relative phases of the two scatterings. 
We find that the emission of photons at about 2 GeV 
is reduced by a factor of 3.5 if the LPM suppression is 
treated in this manner. Most of the enhancement due 
to fragmentation processes visible in Fig. [3] disappears 
in this limit. Also, increasing the pQCD cut-off (p™ m ) 2 
to 1 GeV 2 reduces the number of binary collisions and 
the build-up of cascading activity. The net result is a de- 
crease in the photon yield by a factor of about 1.5 beyond 
Pt = 1 GeV near central rapidities. 
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FIG. 4: Integrated photon multiplicity as a function of N pa , T t. 
Squares denote a calculation with only binary parton-parton 
collisions whereas diamonds show a calculation including frag- 
mentation processes. The scaling with iV*^ t is indicative of 
the photon yield being directly proportional to the number of 
hard collisions. 



functions. We intend to study these effects in a forth- 
coming publication. 

In summary, we have calculated the production of high 
energy photons in collision of gold nuclei at 200 A GeV 
in the Parton Cascade Model. Multiple scattering of par- 
tons is seen to lead to a substantial production of high 
energy photons, which rises further when parton mul- 
tiplication due to final state radiation is included. The 
photon yield is found to scale as and is directly pro- 

portional to the number of hard parton-parton collisions 
in the system, providing valuable information on the pre- 
equilibrium reaction dynamics of the system. Our calcu- 
lations suggest that that at the much higher energies to 
be reached at the Large Hadron Collider, the parton cas- 
cades developing in the nuclear collision will produce an 
even larger "firework" of photons. 

This work was supported in part by DOE grants 
DE-FG02-96ER40945 and DE-AC02-98CH10886 and the 
Natural Sciences and Engineering Research Council of 
Canada. 



How does the photon production scale with the number 
of participant nucleons N palt 7 Figure 0] shows the total 
integrated photon yield (neglecting hadronic and thermal 
contributions) as a function of N paI t. The squares denote 
a calculation without fragmentation processes, whereas 
the diamonds show the full calculation, including frag- 
mentation. For both calculations we observe a scaling of 
the photon yield with N^t (dotted curves), suggestive 
of a scaling with the number of binary parton-parton col- 
lisions 0, [r| rather than participant number. This is a 
valuable result, since it provides a direct connection be- 
tween a measurable quantity (the photon yield) and the 
unobservable number of hard collisions during the time- 
evolution of the reaction. It will be interesting to see in 
future studies whether the production of heavy quarks 
and the suppression of high-px hadrons reflect similar 
scaling properties. 

In our calculation we have used lowest order pQCD 
matrix elements to study the collisions. One could in 
principle use a if -factor to account for higher order ef- 
fects. We have not done so here, as the parton shower 
mechanism for the final state radiations amounts to the 
inclusion of collinear emissions to all orders. At least for 
photon production in a QGP, rates are now available up 
to two loops 0], and to complete leading order [l^. It 
is not clear how these calculations can be incorporated in 
the PCM, unless a (local) temperature is associated with 
the system. Our current calculation neglects medium 
modifications of the matrix elements and fragmentation 
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